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Preface 

In  tills  report,  I  have  attempted  to  present  the  reader 
with  an  easily  understood  report  on  the  evaluation  of  thrust 
coefficient  for  clustered  nozzles* 

This  thesis  represents  a  summary  of  various  concepts 
that  I  have  learned  about  clustering  nozzles.  Although  a 
thesis  is  important  to  the  research  and  scientific  world,  I 
feel  that  my  experiences  and  achievements  exceeded  the  ideas 
contained  'within  these  pages.  During  this  18  month  period, 

I  have  had  the  opportunity  and  pleasure  of  having  a  new  friend. 
As  a  thesis  advisor,  Dr.  W.  C.  Elrod  has  given  of  himself  un¬ 
selfishly.  He  has  always  had  an  ear  to  lend  and  has  offered 
alternative  and  valuable  solutions  to  various  problems  both 
on  the  professional  and  personal  level.  I  have  thoroughly 
enjoyed  working  with  and  having  Dr.  Elrod  as  my  advisor  during 
the  course  of  this  study. 

I  wish  to  acknowledge  my  love  and  appreciation  to  my 
wife,  Kathy.  Throughout  this  study,  I  have  tormented  her  with 
a  sporadic  meal  schedule,  a  never  ending  laundry  pile,  and  a 
mind  plagued  with  engineering  problems.  Yet  she  endured  each 
day  by  handing  me  a  beer  and  telling  me  she  understood.  Also, 
a  special  thanks  to  our  dog  Clyde  who  greeted  me  each  day 
with  a  friendly  wag.  When  things  didn't  go  well  during  the 
course  of  this  study,  as  they  often  didn’t,  Clyde  would  ap¬ 
proach  with  a  ball  to  be  thrown  or  with  a  leash  for  a  walk 
in  the  park.  Mostly,  I  thank  them  both  for  showing  me  in  their 
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own  way  that  John  Lennon  was  right  when  he  said,  "All  you  need 
is  love." 

A  very  special  thanks  to  Dr.  P.  Torvik  for  his  assist¬ 
ance  with  the  visicorder  and  instrumentation.  Professor  Torvik 
made  numerous  and  valuable  suggestions  in  setting  up  and  de¬ 
bugging  the  instrumentation.  I  thank  him  for  taking  a  personal 
interest  in  this  study. 

Many  others  have  contributed  to  this  study  over  the 
past  year.  I  would  like  to  take  the  time  to  thank  a  few  of 
them  now. 

-Dr.  H.  Wright  and  Dr.  M.  Franke,  my  thesis  committee, 
for  their  suggestions  throughout  the  course  of  this 
study. 

-Mr.  Shortt  and  Mr.  Brohas,  from  the  fabrication  shop, 
for  their  handiwork  in  developing  and  building  of  the 
apparatus. 

-Mr.  Baker  and  Mr.  Cannon,  the  lab  technicians,  for 
their  assistance  in  the  assembling  of  the  apparatus. 

-Ms.  K.  Newman,  personal  friend,  for  her  assistance 
with  the  viewgraphs. 

David  V.  Hibson 
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Abstract 

This  is  an  experimental  evaluation  of  the  thrust  per¬ 
formance  of  3  sets  of  clustered,  three-dimensional,  conver- 
gii\g-diverging,  cold  flow,,  supersonic  nozzles*  A  cluster  of 
2,4,  and  6  nozzles  were  designed  and  fabricated.  Each  clus¬ 
ter  assembly  has  the  same  geometry  in  that  their  area  ratio, 
expansion  ratio,  and  total  throat  area  is  the  same.  A  single 
nozzle  with  the  same  geometry  and  total  throat  area  was  used 
to  evaluate  the  creditability  of  the  testing  procedure  and  the 
performance  of  the  3  sets  of  clustered  nozzles.  The  thrust 
performance  of  each  nozzle  cluster  was  evaluated  by  comparison 
of  the  measured  thrust  coefficient  of  the  cluster  to  that  of  a 
single  nozzle.  The  nozzle  with  the  highest  thrust  coefficient 
was  the  cluster  of  2  nozzles.  Its  performance  was  closely 
followed  by  the  cluster  of  4  nozzles.  The  nozzle  with  the 
lowest  thrust  coefficient  was  the  cluster  of  6  nozzles.  The 
results  of  this  study  indicate  that  the  clustering  of  nozzles 
improves  the  thrust  performance. 
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PERFORMANCE  CHARACTERISTICS  OF  CLUSTERED  NOZZLES 


I.  Introduction 

The  concept  of  clustering  arrays  of  rocket  engines  is 
not  new.  In  the  early  60's  nozzles  were  clustered  in  order 
to  capitalize  on  the  altitude  compensating  characteristics  of 
plug  nozzles.  Recent  work  has  been  aimed  at  achieving  high 
vacuum  specific  impulse  by  clustering  existing  rocket  engines 
to  for  very  high  area  ratio  engines.  A  specific  application 
now  being  considered  at  AFRPL,  Edwards  AFB,  California,  is  the 
clustering  of  existing  modules  such  as  the  RL-10  0^/H^  engine. 
No  information  is  presently  available  on  the  performance  of 
clustered  nozzles  of  this  type  at  the  high  pressure  ratios 
characteristic  of  these  rocket  engines  (ref  6). 

The  object  of  this  thesis  wu... 

1.  To  design,  fabricate,  and  calibrate  the  necessary 
apparatus  and  instrumentation  in  order  to  determine 
the  thrust  coefficient  for  single  and  clustered 
nozzles. 

2.  To  experimentally  determine  the  thrust  performance 
and  stability  of  operation  of  3  sets  of  clustered, 
three-dimensional,  converging-diverging,  cold  flow, 
supersonic  nozzles. 

This  is  an  experimental  study  in  model  scale.  In  this 
study,,  a  clustered  nozzle  assembly  is  considered  to  be  2  or 
more  nozzles  that  are  held  together  and  have  equivalent 
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throat  areas.  Determination  of  the  relative  thrust  efficiency 
of  each  nozzle  set  was  accomplished  by  comparison  of  its  thrust 
coefficient  to  that  of  a  single  nozzle  whose  throat  area  was 
the  same  as  that  of  the  total  throat  area  of  the  set. 
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II 


An  objective  of  this  study  was  to  investigate  exper¬ 
imentally  the  thrust  performance  of  clustered  nozzles  in 
model  scale.  The  thrust  performance  of  each  nozzle  cluster 
was  evaluated  by  comparing  it  to  a  standard  single  nozzle 
by  comparing  their  respective  measured  thrust  coefficient. 
The  thrust  coefficient  determines  "the  amplification  of 
tnrust  due  to  the  gas  expansion  in  the  rocket  nozzle  as 
compared  to  the  thrust  that  would  be  exerted  if  the  chamber 
pressure  acted  over  the  throat  area  only"  (ref  5)«  Thrust 
coefficient  may  be  defined  as: 


T 


By  use  of  conservation  momentum,  the  thrust  is  equal 
to  (ref  5) : 

T  —  mV2  ■+■  (p2  “  p^ )  A 2 

"The  thrust  is  composed  of  two  terms.  The  first  term, 
the  momentum  thrust, is  the  product  of  the  gas  (propellant) 
mass  flow  rate  and  the  exhaust  velocity  relative  to  the  noz¬ 
zle  (vehicle).  The  second  term,  the  pressure  thrust,  con¬ 
sists  of  the  product  of  the  cross-sectional  area  of  the  ex¬ 
haust  jet  leaving  the  nozzle  and  the  difference  between  the 
exhaust  pressure  and  the  fluid  pressure.  If  the  exhaust  pre- 
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sure  is  less  than  the  surrounding  fluid  pressure,  the  pressure 
thrust  is  negative"  (ref  5).  This  condition  is  known  as 
overexpansion  and  is  undesirable  due  to  the  loss  in  thrust. 

An  underexpanding  nozzle  discharges  the  gas  at  a  pressure 
greater  than  the  ambient  pressure.  In  this  case,  the  nozzle  exit, 
area  is  too  small.  A  nozzle  with  optimum  expansion  is  when 
the  nozzle  exit  pressure  is  equal  to  the  ambient  pressure.  At 
this  condition,  the  thrust,  and  therefore  the  thrust  coefficient 
also,  is  at  a  maximum. 

The  effect  of  either  overexpansion  or  underexpansion 
is  a  slight  reduction  in  the  exhaust  velocity  and,  therefore, 
a  loss  in  energy.  The  loss  of  thrust  due  to  overexpansion  and 
underexpansion  may  be  determined  from  the  thrust  coefficient. 

This  theory  assumes  that  the  nozzle  is  flowing  full  (i.e. 

P2-  S  0.4p3)  and  that  the  pressure  thrust  is  positive  for 
underexpansion  and  negative  for  overexpansion.  For  under¬ 
expansion  and  for  slight  overexpansion ,  this  simple  theory  is 
in  accurate  agreement  with  measured  results. 

A  sample  calculation  may  be  found  in  Appendix  C.  The 
sample  calculation  gives  rise  to  the  desired  quantities  to  be 
measured.  Knowing  the  desired  quantities,  the  first  objective, 
design,  fabricate,  and  check  the  apparatus,  was  then  achieved. 
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III.  Apparatus 


The  purpose  of  the  apparatus  was  to  provide  a  means 
to  obtain  the  necessary  data  to  determine  the  thrust  coeffi¬ 
cient  for  the  nozzle  clusters.  The  theory  section  indicated 
the  various  quantities  needed  to  be  measured  (i.e.  pressure, 
temperature,  thrust).  From  this,  the  first  objective,  to  de¬ 
sign  and  build  an  apparatus,  was  completed. 

The  experimental  apparatus  consisted  of  a  nozzle  clus¬ 
ter  assembly,  stilling  chamber,  mass  flow  meter,  visicorder, 
and  various  regulating  valves  and  hardware.  Descriptions  of 
eacn  and  its  use  in  this  study  are  given  below.  A  schematic 
is  shown  in  Fig.  1.  Further  details  and  drawings  may  be  found 
in  Appendix  B. 

Nozzle  Cluster  Assembly 

The  nozzles  were  designed  and  then  fabricated  from  aluminum 
stock  (Fig  3)»  The  nozzles  were  designed  using  one-dimensional 
isentropic  relations  (ref  2)  for  a  pressure  ratio  of  £3  =  0.04711 

n 

(M  =  2.64)*  Each  nozzle  inlet  was  designed  using  a  constant 
arc  of  radius  1.837  +0.001  in  (ref  3). 

The  3  sets  of  nozzle  clusters  were  designed  with  the 
same  total  throat  area.  A  cluster  of  2,  4>  and  6  nozzles  were 
assembled.  The  cluster  arrangements  are  shown  in  Fig  4. 

Each  nozzle  was  machined  separately  and  then  fastened  to  an 
aluminum  base  plate.  Clay  was  used  to  seal  the  seamlines  be¬ 
tween  the  nozzles.  For  control  purposes  (see  test  procedure 
section),  a  single  nozzle  with  the  same  geometry  and  total 
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Fig  2  Thrust  Measuring  Setup 
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throat  area  was  made.  Also,  a  spare  nozzle  from  each  cluster 
assembly  was  fabricated.  This  made  a  total  of  7  nozzles,  that 
is  3  cluster  nozzle  assemblies,  a  large  single  nozzle,  and  3 
more  single  nozzles  from  each  cluster. 

Stilling  Chamber 

The  nozzle  base  plate  fastened  directly  to  the  stilling 
chamber.  The  chamber  provided  a  reservoir  of  compressed  air. 
Within  the  stilling  chamber,  a  basket-type  diffuser  (Fig  B-l) 
disseminated  the  incoming  gas  in  a  radial  direction  to  the 
stilling  chamber.  The  nozzles  were  tested  at  a  chamber . pressure 
of  150,  200,  250,  300,  and  350  psig. 

The  stilling  chamber  pivoted  via  its  3  ft  long  feed 
line  as  shown  in  Fig  2.  At  the  pivot  point,  the  inner  slotted 
tube  diffused  the  gas  in  a  radial  direction  before  the  gas 
entered  the  feed  line.  Bearings  and  0-rings  were  used  to  re¬ 
duce  friction  and  provide  a  proper  seal  (Fig  B-2).  Strain 
gages  mounted  on  an  aluminum  beam  were  used  to  measure  the 
thrust  created  by  the  nozzle  assembly  (Fig  2).  A  pressure 
transducer  and  a  thermocouple  measured  the  pressure  and  tem¬ 
perature  of  the  gas  respectively. 

Mass  Flow  Meter 

A  thin-plate  square  edge  orifice  meter,  built  and  in¬ 
stalled  according  to  ASME  standards,  was  located  just  prior 
to  the  stilling  chamber  assembly  (Fig  1).  The  mass  flow  meter 
was  designed  using  the  energy  equation  (ref  1)  for  a  diameter 
ratio  of  =  i  (Fig  B-3)»  Two  pressure  transducers  and  a 
thermocouple  measured  the  pressure  and  temperature  of  the  gas. 
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Visicorder 


A  Honeywell  visicorder,  model  1508,  was  used  to  re¬ 
cord  all  data.  The  voltage  outputs  from  the  various  pressure 
tranducers,  thermocouples,  and  thrust  measuring  device  were 
each  wired  independently  to  one  of  the  8  channels  to  the  vi¬ 
sicorder.  Of  the  8  channels  used,  the  thrust  data  was  record 
ed  on  1  of  3  channels.  The  thrust  measuring  channel  was 
selected  according  to  the  nozzle  being  tested  (see  test  pro¬ 
cedure). 

Regulating  Valves  and  Hardware 

A  Grove  regulating  valve  was  used  to  adjust  the  pres¬ 
sure  above  the  diaphram  in  the  dome  regulator.  The  dome  regu 
lator  then  adjusted  the  air  pressure  from  the  supply  line 
pressure  to  the  desired  operating  pressure.  Air  was  supplied 
from  a  nearby  trailer  unit. 

After  the  apparatus  was  completed,  the  instrumentation 
was  calibrated  and  test  procedures  were  established  as  des¬ 
cribed  next. 


IV.  Test  Procedure 


V/ith  the  experimental  apparatus  assembled,  the  last 
part  of  objective' 1,  to  calibrate  the  instrumentation  and 
to  provide  a  test  procedure  for  the  apparatus,  was  completed. 
Nozzles  1  thru  4  were  used  as  a  standard  to  calibrate  and 
de-bug  the  system.  A  description  of  how  the  mass  flow  meter, 
thrust  measuring  device,  pressure  transducers  and  thermo¬ 
couples,  and  visicorder  were  calibrated  are  given  below. 

Mass  Flow  Meter 

The  mass  flow  meter  was  designed  and  installed  according 
to  ASME  standards  (ref  1).  No  direct  calibration  was  required. 
Various  equations  and  charts  (see  ref  1  and  appendix  C)  were 
used  to  determine  the  mass  flow  rate.  As  a  check,  isentropic 
conditions  were  assumed  at  the  nozzle  throat  to  determine  a 
second  calculation  of  mass  flow  rate.  The  2  mass  flow  rates 
were  used  in  the  results  section. 

Thrust  Measuring  Device 

The  thrust  measuring  device  was  an  aluminum  beam  with 
mounted  strain  gages  as  shown  in  Fig  2.  The  thrust  is 
measured  along  the  centerline  of  the  stilling  chamber  and  the 
nozzle.  This  design  restricted  the  movement  of  the  stilling 
chamber  to  5/32  in  and  thereby  reduced  unnecessary  friction 
from  the  O-ring  assembly  at  the  pivot  point.  To  further 
reduce  this  effect,  the  stilling  chamber  was  vibrated  during 


the  data  collecting  period.  This  method  reduced  friction 
effects  to  a  minimum. 

A  similar  problem  was  encountered  during  the  calibra¬ 
tion  of  the  thrust  measuring  device.  . Weights  were  suspended  ■ 
from  the  device  via  a  basket,  cable  and  pulley  assembly 
(Fig  2).  Again  the  stilling  chamber  was  vibrated  to  reduce 
frictional  effects  from  the  pulley  assembly. 

This  study  does  not  concern  itself  with  the  exact  ef¬ 
fects  of  the  friction  incurred  from  the  calibration  and  test 
procedures  on  the  data.  Nozzle  1  was,  tested  to  indicate  if 
the  hysterisis  effects  were  indeed  negligible  (Fig  A-ll).  A 
correlation  factor  was  found  for  the  thrust  data  taken  and 
indicated  that  the  data  was  linear.  Knowing  that  the  mea¬ 
sured  thrust  data  obtained  is  linear,  hysterisis  from  fric¬ 
tional  effects  was  assumed  negligible  in  the  testing  procedure. 
Pressure  Transducers  and  Thermocouples 

Pressure  transducers  were  calibrated  by  use  of  a  dead 
weight  tester.  Two  iron  vs  constantan  thermocouples  were 
used  to  measure  the  gas  temperature.  They  were  calibrated  by 
inserting  the  thermocouple  bead  into  a  water  bath  with  a 
thermometer,  heating  it,  and  recording  the  visicorder  read¬ 
ing  and  the  thermometer  reading. 

ViSLSSXler 

All  measuring  devices  fed  their  voltage  outputs  directly 
to  the  visicorder.  The  visicorder  recorded  all  inputs  on  a 
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light  sensitive  graph  paper.  All  data  was  recorded  with¬ 
in  the  20  cm  width  of  paper.  In  order  to  provide  the  proper 
resolution,  the  thrust  measuring  device  was  required  to 
record  on  3  separate  channels.  One  channel  for  thrust  0  - 
35  lt>£,  another  for  25“85  lb^,  and  the  third  for  40  - 
150  lb^  provided  the  data  to  be  resolved  to  within  2.0% 
of  its  recorded  value. 

All  inputs  to  the  visicorder,  with  the  exception  of 
the  thermocouples,  used  a  variable  cesistor  as  indicated  in 
Fig  3-4.  The  variable  resistor  provided  2  useful  functions. 
First,  ±\  permitted  the  entire  recording  space  to  be  used. 

An  example  of  this  is  the  pressure  transducer  for  the  stilling 
chamber.  Tests  were  made  between  1 50  -  350  psig.  There¬ 
fore,  the  1  cm  mark  was  the  1 50  psig  point  and  the  19  cm 
mark  was  the  350  psig  point.  This  permitted  the  best  at¬ 
tainable  resolution  in  the  recorded  data  and  also  kept  the 
galvanometer  well  within  its  linear  range.  Second,  the  var¬ 
iable  resister  permitted  an  easy  method  to  calibrate  the 
readings  from  day  to  day.  An  example  of  this  would  be  to 
place  a  15  lbf  on  the  thrust  measuring  device.  The  variable 
resistor  could  then  be  adjusted,  if  necessary,  to  locate  the 
light  beam  at  the  proper  point.  In  addition,  the  variable 
resistor  did  not  change  the  sensitivity  of  the  galvanometer. 
Procedure 

A  nozzle  assembly  was  selected  and  fastened  to  the 
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rear  of  the  stilling  chamber  (Fig  2).  The  Grove  regulator 
was  then  adjusted  to  increment  the  gas  pressure  in  the  stil¬ 
ling  chamber  from  150  -  350  psig  and  back  to  150  psig  in  50 
psig  increments.  The  apparatus  was  vibrated  to  reduce 
frictional  effects.  This  method  produced  data  with  a  high 
degree  of  repeatability  as  the  results  section  indicates 
next. 
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V.  Results  and  Discussion 

The  apparatus  was  designed,  fabricated,-  and  checked 
out  with  the  testing  of  nozzle  1.  The  single  nozzle  thrust 
performance  was  compared  to  theory  (Fig  A-ll)  in  order  to 
establish  a  base  line.  A  statistical  package  within  the 
AFIT  computer  system  was  used  to  fi c  the  best  curve  to  the 
measured  and  theoretical  data.  A  linear  curve  was  fitted 
to  the  measured  thrust  data  points.  A  correlation  factor 
of  0.9996  indicated  that  the  line  fitted  to  these  points  was 
nearly  an  ideal  fit.  Knowing  that  the  measured  data  points 
were  linear,  it  was  assumed  that  any  hysterisis  effects 
encountered  in  the  calibration  procedures  would  have  little 
effect  in  the  data  collecting  procedures.  From  the  thrust 
measurements,  the  thrust  coefficient  was  determined  (Fig  A-4). 
Various  statistical  packages  within  the  computer  system  were 
used  to  fit  polynomial  curves  to  the  data.  Nozzle  1  estab¬ 
lished  a  baseline  upon  which  the  3  nozzle  cluster  sets  and 
their  respective  standards  were  evaluated. 

A  typical  graph  of  the  thrust  coefficient  vs  pressure 
ratio  is  shown  in  Fig  5.  Graphs  of  all  measurements  may  be 
found  mn  Appendix  A.  A  sample  calculation  with  an  error 
analysis  may  be  found  in  Appendix  C.  Table  I  shows  the  order 
of  best  performing  nozzles  at  optimum  expansion,  Pr=  21.13. 

The  equations  derived  by  the  computer  may  be  found  at  the 
end  of  Appendix  A.  The  theory  curve  in  Fig  3»  and  A-l  thru 
A-3  was  found  by  using  the  theoretical  data  points  for  all 
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NOZZLE  2-2  *  *  - 

NOZZLE  3-4  i  $  - 

NOZZLE  4-6  »  §  - 

NOZZLE  1  «  «• - 


FIG  5  i  CLUSTERED  AND  SINGLE  NOZZLES 


the  nozzles  in  that  particular  graph.  This  resulted  in 
slightly  different  equation  (due  to  0.7%  difference  in  At) 
as  shown  in  Table  A-2. 

Comparison  of  Cluster  Sets  with  Single  Nozzle 

A  comparison  of  the  3  cluster  nozzle  sets  with  a  sin¬ 
gle  nozzle  of  equivalent  throat  area  is  shown  in  Fig  5* 

From  the  calculations  in  appendix  C,  the  optimum  expansion 
point  is  at  Pr=  21,13.  A  pressure  ratio  less  than  21.13 
operates  the  nozzle  at  an  overexpanded  condition  while 
pressure  ratios  greater  than  21.13  operate  the  nozzle  at 
an  underexpanded  condition.  Nozzle  2-2  has  a  1.6%  increase 
in  thrust  compared  to  the  single  nozzle  at  P  =  21.13. 

Similarly,  nozzle  3-4  has  a  1.0%  increase  in  thrust.  Nozzle 
4-6  v/as  slightly  less  than  the  single  nozzle  by  -0.2%.  Over¬ 
all,  for  all  pressure  ratios,  nozzle  2-2  was  the  best  performer. 
Comparison  of  Nozzle  2-2  with  Nozzle  2 

Nozzle  2-2  is  compared  with  nozzle  2  in  Fig  A-l.  At 
a  Pr=  21.13,  nozzle  2-2  has  a  2.2%  increase  in  thrust  over 
nozzle  2.  For  all  pressure  ratios,  nozzle  2-2  out  performs 
nozzle  2. 

Comparison  of  Nozzle  3-4  with  Nozzle  3 

The  results  of  nozzle  3-4  and  nozzle  3  are  shown  in 
Fig  A-2,  Nozzle  3-4  is  the  best  performer.  At  Pr=  21.13* 
nozzle  3-4  has  a  0.7%  increase  in  thrust  over  nozzle  3. 
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Comparison  of  Nozzle  4-6  with  Nozzle  U 

A  graph  of  the  results  for  nozzle  4-6  and  nozzle  4  is 
shown  in  Fig  ,A-3.  Here,  nozzle  4-6  out  performs  nozzle  4 
until  a  pressure  ratio  of  about  15  is  reached,  and  then  the 
single  nozzle  4  performs  better  than  nozzle  4-6.  At  optimum 
expansion,  nozzle  4  has  a  0.7%  increase  in  thrust  over  nozzle 
4-6. 

Nozzle  1 

The  results  for  nozzle  1  are  plotted  in  Fig  A-4.  At 
P  =  21.13,  it  is  within  3*3%  of  theoretical  results.  All 

tL 

nozzles  are  tabulated  in  Table  I. 


Table  I.  Performance  Ratings  at  Optimum  Condition 


Nozzle 

Cf 

t 

Cf 

m 

A% 

2-2 

1.38 

1.36 

1.5 

3-4 

1.38 

1.35 

2.4 

4 

1.39 

1.35 

2.9 

3 

1.39 

1.35 

5.1 

1 

1.39 

1.34 

3.3 

4-6 

1.38 

1.34 

3.4 

2 

1.39 

1.34 

3.5 

A%  = 


Cf  -  c 
xt  xm 


*  100  % 


m 


18 


-U— A. 


i  *  "7  ?  3"  ,*  ,•  ^rrrrr  -  ^ ^ 

4?  ' d'--.:  ;•  ’It  ’a  i  uV- 


VI.  Conclusions 

The  apparatus  and  instrumentation  was  designed, 
fabricated  and  check  out  by  the  author.  From  the  results 
section  and  the  error  analysis  in  Appendix  C,  the  data 
obtained  from  this  study  has  a  high  degree  of  repeatability. 
The  error  analysis  indicated  that  the  data  points  should 
have  a  spread  of  +  2,0%,  The  actual  data  points  are  con¬ 
tained  well  within  these  boundaries.  From  this,  the  fol¬ 
lowing  conclusion  may  be  drawn: 

1.  The  apparatus  is  suitable  for  determining  the 
thrust  coefficient  for  single  and  clustered  nozzles. 

A  single  nozzle  was  tested  and  evaluated  in  order  to 
create  a  baseline.  The  3  sets  of  clustered  nozzles  were  then 
compared  to  the  single  nozzle.  Following  this,  the  3  nozzle 
sets  and  their  respective  standard  nozzles  were  evaluated. 
After  analyzing  the  results,  the  following  conclusion  may 
be  drawn: 

2.  The  thrust  performance  does  vary  with  each  of  the 
different  nozzle  configurations.  This  study  suggests  that 
there  is  a  relationship  between  the  number  and  arrangement 
of  the  nozzle  exhausts  and  the  free  stream  flow  such  that 
their  interaction  influences  the  pressure  thrust  term  in  the 
thrust  equation.  This  study  indicates  that  clustering  of 
nozzles  improves  thrust  performance. 
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The  recommendations  arising  out  of  the  course  of  this 
study  are: 

1.  Further  study  of  clustered  nozzles  should  be 
performed  to  examine  the  effect  of  shrouding  between  the 
nozzles  to  improve  the  performance  of  the  4  and  6  nozzle 
cluster. 

2.  Additional  study  of  clustered  nozzles  of  non- 
convensional  shape  (i.e.  rectangular,  oval,  etc.)  and  in 
various  configurations  could  be  completed  to  provide  a 
data  base  for  future  design  work. 

3.  Continue  this  study  using  higher  pressure  ratios. 
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Table  A-i  Nozzle  Equations  Derived  by  Computer  for  Fig  A-1 
' hru  A- 17  and  nut  into  Tex 


Nozzle 

y 

a(2)PIT 

+  a  ( 1 )  P 

P  +  b 

2-2 

% 

-0.0006 

0.0329 

0.9662 

-0.0004 

0.0236 

1 .0260 

Tt 

6.45 

-12.86 

T 

m 

6.33 

-12.93 

3-4 

V 

-0.0006 

0.0331 

0.9633 

-0.0007 

0.0366 

C.9103 

Tt 

6.44 

-12.72 

T 

m 

6.28 

-12.24 

4-6 

Cf 

2t 

-0.0006 

0.0323 

0.9738 

Cf 

m 

-0,0006 

0.0319 

0.9346 

Tt 

6.45 

-12.71 

Tm 

6.23 

-12.18 

1 

% 

-0.0007 

0.0371 

0.9323 

-0.0006 

0.0318 

O.929I 

-T- 


Table  A-2  Averaged  Theory  Equations  Derived  by  Computer  and 


a(2)PIT  + 

a(1 )PR 

+  b 

-0v,000<7 

0.0351 

0.9466 

-0.0005 

0.0293 

0.99 28 

-0,0006 

0.0318 

0.9754 

-0.0006 

0.0335 

0.9616 

I k«»'  ir?i**  i  - 


voltage 

supply 


Fig  B-4  Variable  Resistor  Setup 
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Sample  Calculation 


The  thrust  coefficient  determines  the  amplification 
of  thrust  due  to  the  gas  expansion  in  the  rocket  nozzle  as 
compared  to  the  thrust  that  would  be  exerted  if  the  chamber 
pressure  acted  over  the  throat  area  only  (ref  5) •  The  basic 
definition  of  thrust  coefficient  is: 


where:  T  =  thrust  (lbf) 

A^  =  throat  area  (in^) 

Pi  =  chamber  pressure  (psia) 

Thrust  coefficient  is  a  non-dimensional  term  that  is 
used  in  this  study  to  compare  the  thrust  performance  among 
the  7  nozzles  tested.  The  measured  thrust  coefficient  was 
found  by  measuring  the  thrust,  throat  area  and  chamber  pres¬ 
sure  directly.  To  determine  the  theoretical  thrust  coefficient, 
more  extensive  calculations  were  carried  out.  To  obtain  the 
the  theoretical  thrust,  Sutton  (ref  5)  used  conservation  of 
momentum  to  find: 


(  ~  )  A ^ 


49 


where:  T  =  thrust  (lbf) 

m  =  mass  flow  rate  (lb  /s) 

m 

V2  =  exit  velocity  of  gas  (fps) 
ft' lb  /s2 

So  ■  "lb”1 

P2  =  pressure  at  exit  plane  of  nozzle  (psia) 

=  ambient  pressure  (psia) 

p 

A2  =  nozzle  exit  area  (in  ) 

A2  can  be  measured  directly  from  the  nozzle  itself  and 
is  the  barometric  pressure  during  the  testing  period.  Vp 
and  p2  are  found  using  isentropic  relations  from  ref  2.  The 
mass  flow  rate  (m)  can  be  determined  from  two  methods. 

Method  I  ;  .Mass  Flow  Me.t.e.r. 

The  following  equations  were  obtained  from  ref  1  (p58) 
and  were  derived  from  the  energy  equation. 

D  =  1.473”  8  =  0.737" 

ft  =  d/D  =  0.500 

K  =  0.5993  +  0.007/D  +  (0.564  +  0.076 *TB) 

V 

K  =  0.6308 

A  =  d(830  -  5000^  +  9000^ 2  -  420(jtf3  +  530/Vd’  ) 

A  =  362.59 

K  =  Ko(d*106/(d*106  +  1 5A) ) 

U  c 

KQ  =  0.6262 

50 


K  =  C/V  (1  -  /sF)  =  K0(1  +  VRd) 

Rd  =  48m/('rrd/i) 

where:  m  =  1.37  lb  /s  (assumed) 

JJ=  1.19923  105  lbffl/ft  s 

Rd  =  2.3  106 

K  =  0.6262 

then: 

a  =  <r/4*d2*K*  ^(2gc/oa(pa  -  Pb)>' 

A  =  V(z  K  Ta) 

where:  z  =  compressibility  factor 

=  0.99  @  Pa  =  330  psia  Ta  =70F 

R  =  gas  constant 

=  55.33  lbf  ft/ (R  lbm) 

substituting: 

a  =  0.2895  V  (Pa/Ta)(pa  "  pb}’  lbm/s 

where:  p  =  pressure  before  orifice  (psia) 
SI 

pd  =  pressure  after  orifice  (psia) 


=  temperature  of  gas  (R) 


4  *  AAtvt 


for  choked  flow  M  =  1 

ref  .2  Tt  =  0.8333  T1 

Pt  =  0.528282  p1 


V,  =  M 


t  ttt 


where:  =  1.0 

at  »  o*RTt 


A  =  V(KV 


where:  R  =  53.33  lof  ft/(R  lbm) 


substituting: 

m  =  0.522  (Pl  A^/T^ )  lt>m/s 

where:  p^  =  chamber  pressure  (psia) 
=  chamber  temperature  (R) 

p 

A^  =  throat  area  (in  ) 


Continuing  with  the  thrust  equation: 


A2  is  measured  directly  from  the  nozzle  exit. 

A2  =  'n‘/4*d2  in2 

The  exit  areas  were  as  follows: 

1  0.8882759  in2 

2  0.4624008  in2  2-2  0.8835211  in2 

3  0.2206246  in2  3-4  0.8824986  in2 

4  0.1520122  in2  4-6  0.8835729  in2 
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The  thrust  can  now  be  calculated  for  an  assumed  and 
Tj.  It  is  assumed  that  the  temperature  measurement  is  reasom- 
ably  correct  and  can  be  used  for  the  theoretical  calcula¬ 
tion.  Temperature  of  the  gas  was  measured  in  degree  F.  V/hen 
converted  to  absolute  temperature  and  the  square  root  is 
taken,  a  temperature  uncertainty  of  +  I.OF  has  little  effect 
in  the  theoretical  calculation. 

Also  it  should  be  noted  that  optimum  expansion  takes 
place  when  p£  =  p^.  This  occurs  when: 

=  1/0.047329 
*3 

or  Pp  =  21.129 

The  theoretical  and  measured  thrust  coefficients  can 
now  be  determined: 
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firm  i&aijsls 

In  measuring  the  data,  the  visicorder  could  give 
the  following  resolutions  in  its  recordings: 

P1  =  150  +  1.25  psig 
T  =  10  +  0.175  lb^  for  nozzles  3  &  4 
=  35  +  0,28  lbf  for  nozzle  2 
=  60  +  0,59  lbf  for  nozzles  1,  2-2,  3-4,  &  4-6 

In  measuring  the  throat  areas,  the  diameters  could 
be  measured  to  within  +  0.0001  in. 

=  0.0507  +  0.00016  in2  for  nozzle  4 
=  0.1541  +  0.00028  in2  for  nozzle  2 

M 

=  0.2961  +  0.00037  in2  for  nozzle  1 

r  _ T _ 

Cf  -  At?1 

=  1 

d  T  AtP1 


Given: 


then: 
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fif' 


•.Vj'ft.r 


-  L  (.  _i)2 

"  At  pi 


^Cf  iL  f  J-\2 

■  pT  ('  V 


then: 


(5t*  <^T)2  +  (|~  ^  )2  +  (jlf 


_  +  r  ( _ Q.tiz5„,.)2  +  ( _ ici 

-  i  |_m  50*0. 0507 }  0.050*; 


10*1.25  s2 


7*150 


+  ( 10*0.00016  n2 

150*0. 05072  J 


7a 


=  0.026 


Therefore:  Cf  =  Cf  +  0.026 


for  nozzle  3  &  4 


Similarly:  +  0.018  for  nozzle  2 


Cf  =  Of  +  0.017 


for  nozzle  1,  2-2,  3-4, 

&  4-6 


From  this  error  analysis,  a  region  within  each  data 
point  may  be  drawn  to  indicate  a  region  of  uncertainty.  All 
of  the  data  points  for  each  graph  were  used  in  order  to  fit 
the  best  curve.  In  comparing  the  data  point's  region  of  un¬ 
certainty  to  its  vicinity  to  the  curve,  the  data  points  indi- 

a  high  degree  of  repeatability  from  which  trends  may  be  shown 

in  the  nozzle’s  thrust  performance. 
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